On the jets, kinks, and spheromaks formed by a planar magnetized coaxial gun 
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Measurements of the various plasma configurations produced by a planar magnetized coaxial 
gun provide insight into the magnetic topology evolution resulting from magnetic helicity injection. 
Important features of the experiments are a very simple coaxial gun design so that all observed 
geometrical complexity is due to the intrinsic physical dynamics rather than the source shape and 
use of a fast multiple-frame digital camera which provides direct imaging of topologically complex 
shapes and dynamics. Three key experimental findings were obtained: (1) formation of an axial 
coUimated jet [Hsu and Bellan, Mon. Not. R. Astron. Soc. 334, 257 (2002)] that is consistent with 
a magnetohydrodynamic description of astrophysical jets, (2) identification of the kink instability 
when this jet satisfies the Kruskal-Shafranov limit, and (3) the nonlinear properties of the kink 
instability providing a conversion of toroidal to poloidal flux as required for spheromak formation 
by a coaxial magnetized source [Hsu and Bellan, Phys. Rev. Lett. 90, 215002 (2003)]. A new 
interpretation is proposed for how the n = 1 central column instability provides flux amplification 
during spheromak formation and sustainment, and it is shown that jet coUimation can occur within 
one rotation of the background poloidal field. 

PACS numbers: 52.55.Ip,52.72.+v,52.35.Py,95.30.Qd 



INTRODUCTION 



There has long been a mutually beneficial exchange of 
ideas between laboratory plasma physics and plasma as- 
trophysics. An example is the proposal by Lundquis^ 
that cosmical magnetic fields are force-free^i satisfying 
V X B = AB. Woltjer showed^^ that these force-free 
states are intimately related to the concept of magnetic 
helicity conservation and that the force-free state results 
from a variational principle in which the magnetic en- 
ergy of a system is minimized subject to the constraint 
of constant magnetic helicity, which is a measure of mag- 
netic flux twist and linkage."* This idea was then used 
to explain the spontaneous generation of reversed mag- 
netic fields in the reversed field pinch^ and the tendency 
of coaxial gun produced plasmas to reach the spheromak 
state. ^ Although Taylor relaxatiorJ as this process is now 
known, is not a complete and rigorous theory, it was very 
successful, indicating that the same fundamental physical 
principles govern both laboratory and astrophysical plas- 
mas. The Taylor argument showed that the end result of 
complex dynamics is a relatively simple relaxed state, but 
the argument sidesteps the issue of characterizing the de- 
tailed dynamical processes leading to this state. This pa- 
per describes recent experiments providing new insights 
into the detailed dynamical processes of magnetic helic- 
ity injection and plasma relaxation. In particular, the 
experiments shed light on the non-axisymmetric dynam- 
ics of spheromak formation, a subject with a long his- 
tory beginning with the work of AlfvenS^ and Lindberg,^ 
These experiments also show how magnetic helicity in- 
jection can lead to the formation and coUimation of 
magnetically-driven jets in astrophysics. ■'•SiiLiS 

An important feature of this experiment is the use of 
the simplest possible setup for allowing studies of the de- 
tails and three-dimensional (3D) dynamics of magnetic 



helicity injection and spheromak formation. Key results 
include (1) observation of distinct plasma regimes de- 
pending on the parameter Agun (ratio of peak gun cur- 
rent to flux) , similar to those previously foundM using a 
conventional cylindrical coaxial gun, (2) formation and 
axial expansion of a central plasma column from merg- 
ing of plasma-filled poloidal fiux tubesi^ in a manner 
consistent with magnetohydrodynamic (MHD) models of 
magnetically-driven astrophysical jets^SiiSiiSiii (3) clear 
identification of a kink instability of this central column 
and a kink onset threshold that is quantitatively consis- 
tent with the Kruskal-Shafranov limit for an ideal MHD 
kink,^**'*^ and (4) identification of the nonlinear limit of 
the kink instability as a poloidal flux amplification mech- 
anism leading to spheromak formationsiS, 

This paper is organized as follows. Sectionimdescribes 
the experimental setup. Section IlIII presents the experi- 
mental results. Section ITVI provides a detailed discussion 
of the results in the context of both spheromak formation 
and astrophysical jets. A summary is given in Sec. 



II. EXPERIMENTAL SETUP 

In order to facilitate interpretation of results, the coax- 
ial gun design was made to be as geometrically simple 
as possible: a disk surrounded by a co-planar annulus 
(Fig. 1) and a bias magnetic field coil located just be- 
hind the disk (not visible in Fig. 1). This experimental 
setup provided the additional advantage that the entire 
plasma formation region is in plain view so that all steps 
of the topological evolution coifld be observed. Note that 
this setup has essentially the same magnetic topology as 
a star-disk system associated with magnetically-driven 
astrophysical jets (see Fig. 2). The various components 
of the experimental setup are described in detail below. 
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A. Vacuum chamber 

The cylindrical vacuum chamber, shown in Fig. 3 (top), 
is about 2 m long, 1.4 m in diameter, and constructed 
using 6 mm thick electro-polished stainless steel. The 
chamber has large windows on the side that provide a 
good view of the entire plasma. A key point is that the 
chamber size is large compared to both the coaxial gun 
source and the plasma, so that plasma-wall interactions 
are minimal during the time of interest for the results pre- 
sented here. The chamber is evacuated to 1.5 x 10~^ T 
with a cryogenic pump, and a sorption pump is used for 
initial pump-down from atmosphere to 100 mT so that 
the the entire pumping system is oil-free. The oil-free 
property and the minimal plasma-wall interactions re- 
sulted in plasmas that did not require a discharge clean- 
ing protocol in order to provide clean, reproducible re- 
sults. 



B. Planar magnetized coaxial gun 

The planar magnetized coaxial gun^S, shown in Figs. 1 
and 3 (bottom), is mounted on one end-dome of the vac- 
uum chamber, as shown in Fig. 3(top). The gun consists 
of (1) an inner 20.3 cm diameter copper disk electrode, 
(2) a co-planar, coaxial outer 50.8 cm diameter copper 
annulus electrode, and (3) a magnetic field coil located 
just behind the inner electrode. The copper disk inner 
electrode is bolted onto the vacuum-face of a re-entrant 
7.5 in. diameter stainless steel half nipple mounted on 
one end of an 8 in. diameter ceramic break which insu- 
lates the inner electrode from ground. The other end of 
the ceramic break is mounted on a 10 in. flange centered 
on the vacuum chamber end-dome. Thus, the inner elec- 
trode is electrically floating and is connected to a capaci- 
tor bank via coaxial cables connected across the ceramic 
insulating break. The outer annulus electrode is mounted 
via copper brackets directly to the inner surface of the 
vacuum chamber end-dome and therefore is electrically 
tied to chamber ground. A 6 mm vacuum gap separates 
the disk and outer annulus. The experiments are char- 
acterized using the cylindrical coordinate system shown 
in Fig. 3(bottom); the </> and (R-Z) directions will be re- 
ferred to as toroidal and poloidal, respectively. The disk 
and annulus each have eight 6 mm diameter gas injection 
holes spaced equally in cj). The gas injection holes on the 
disk and annulus are at i? = 4.8 and 17.8 cm, respec- 
tively. Details of gas injection are given in Sec. Ill Dl The 
external bias magnetic field coil is located just behind 
the central disk but, because of the re-entrant geome- 
try, is at atmosphere. The coil consists of 110 turns of 
#11 insulated square magnetic wire wound on a phenolic 
form (13.0 cm inner and 18.4 cm outer diameters) and 
has an inductance of 2.3 mH and resistance of 0.39 fl. 
Contours of constant poloidal flux ip generated by the 
coil are shown in Fig. 3 (bottom). 



C. Power systems 

The coaxial gun is powered by an ignitron-switched 
120 fjF capacitor bank rated up to 14 kV. Eight 2 m long 
Belden YK-198 low inductance coaxial cables in parallel 
carry power from the capacitor bank to the gun elec- 
trodes. Routine operation at bank voltages of 5 kV yields 
peak plasma currents /gun ~ 80 kA and gun voltages 
V^un ~ 2.5 kV (after breakdown), or peak gun power on 
the order of 200 MW. Typical /gun and Vgun traces are 
shown in Fig. 4. 

The external magnetic field coil is powered by a 
1.4 mF, 450 V capacitor bank. The half-period of the 
current trace for the external coil is 10 ms and is effec- 
tively constant on the 10 /iS time-scale of the experiment. 
The bias poloidal magnetic flux i/'gun (intercepting the in- 
ner gun electrode) created by the external coil is typically 
1-2 mWb. 



D. Gas injection 

Fast gas puffs are required in these experiments be- 
cause local gas pressures on the order of 100 mT are 
required for breakdown. A slow or steady gas flU is in- 
tolerable because the entire vacuum chamber would rise 
to 100 mT which would result in a weakly ionized cold 
plasma and also an overloading of the cryopump. Fast 
gas valves provide a transient, highly localized cloud of 
high pressure gas in front of the gun electrodes. The 
gas valves utilize a pulsed current in a thin pancake coil 
to induce image currents in an adjacent aluminum diskSS 
but also employ a plenum for measured, reproducible gas 
output. The disk is thus repelled from the coil, creating 
a transient opening for gas flow from the plenum into 
the vacuum chamber. Restoring force on the disk is pro- 
vided by a combination of the high pressure gas in the 
gas lines and a metal spring. Calibrations indicate that 
each pulse injects lO^'^-lO^^ hydrogen molecules, depend- 
ing on the applied voltage. Optimum timing of gas valve 
flring is determined empirically by adjusting the valve fir- 
ing time to minimize the delay between capacitor bank 
trigger and gas breakdown (typically a few /is). One 
valve is used for the eight inner electrode gas injection 
holes, and two valves are used for the eight outer elec- 
trode holes. Typical gas valve firing voltage is 600 V, with 
all three valves injecting on the order of 10^^ hydrogen 
molecules. It was found that the amount of gas injected 
per pulse slowly increased throughout the dajiSi due to 
heating in the electrolytic capacitors^ powering the gas 
valve. A new power supply using metal film polypropy- 
lene capacitors has been constructed recently to eliminate 
this problem>2^ 
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E. Diagnostics 

The images presented in this paper were taken with 
muhiple-frame charge-coupled device (CCD) cameras. 
Most of the images were taken with a Cooke Co. HSFC- 
PRO (Figs. 5, 8, 12, and 14). The specifications for this 
camera are as follows: 4 or 8 frames per shot, 1280 x 1024 
pixels per frame, 12 bits per pixel, a 10 ns minimum ex- 
posure time, and a minimum inter-frame time of 5 ns 
(250 ns) in 4 (8) shot mode. Other images were taken 
with a DRS Hadland Imacon 200 (Fig. 11), which has 8 
or 16 frames per shot, 1200 x 980 pixels per frame, 10 
bits per pixel, a minimum 10 ns exposure time, and a 
minimum of 5 ns (250 ns) inter- frame time in 8 (16) shot 
mode. The camera is placed on a tripod and positioned 
in front of the right-most vacuum chamber window shown 
in Fig. 3(top). 

Magnetic field data are taken with a radial array of 
small pickup coils mounted on a stainless steel shaft, 
which can be scanned in Z and rotated in (\). The array 
contains twenty groups of three coils arranged so that 
all three components of B are measured from i? = to 
i? = 40 cm with 2 cm resolution. The coils are commer- 
cial chip inductors (Coil Craft 1008CS-472XGBB) and 
have calibrated turns x area N A « 1.3 cm^. The coils 
are placed into precision-machined slots of a long thin 
strip of Delrin, which in turn is slid down a thin-wall 
stainless steel tube covered by an alumina tube to pre- 
vent metal contact with the plasma (see Fig. 1). The 
alumina has an outer diameter of 8.4 mm. The frequency 
response, limited by the skin effect of the stainless steel 
tube, is good up to 1 MHz. Further details of probe de- 
sign, construction, calibration, as well as probe improve- 
ments since this work are reported elsewhere»2^ 

A Rogowski coil placed around the ceramic break, 
which is connected to the inner electrode, measures total 
/gun- The Rogowski coil was hand- wound with a cali- 
brated N A = 9.913 x 10~^ m^, and the signal is inte- 
grated via a passive integrator with RC = 8 ms. An at- 
tenuating Tektronix P6015 high voltage probe measures 

F. Control and data acquisition 

Triggering of the various power supplies (bias coil, 
gas valves, and coaxial gun), the data acquisition sys- 
tem, and the camera is provided by a programmable 12- 
channel sequencer with fiber optic outputs. A typical 
trigger sequence is as follows: (1) the bias coil is fired 
at t — —10 ms (due to its w 10 ms rise-time), (2) the 
gas valves are fired at t = —1.5 ms (to accommodate the 
« 1.5 ms travel time of neutral gas through the gas lines 
into the plasma formation region), and (3) the main ca- 
pacitor bank ignitrons are triggered at i = with gas 
breakdown typically occuring at i « 4 /is. 

The digital data acquisition system is a 64-channel 
VME (Versa Module Europa) system from Struck Inno- 



vative Systems. Each channel has a sampling rate of 
105 MHz, a 12-bit resolution, 50 f7 input impedance, 
32 kB memory, and an input range of -512 to 512 mV 
with continuously adjustable offset. All magnetic probe 
and gun diagnostic signals are digitized on this system. 
The digital acquisition system is controlled and the data 
is displayed and analyzed using Interactive Data Lan- 
guage (IDL) routines. The magnetic probe array signals 
are integrated numerically. 

G. Plasma parameters 

The plasmas produced by the coaxial gun have the 
following nominal global parameters: n ^ 10^'^ cm~^^ 
^ Ti 5-15 eV, and B « 0.2-1 kC. The global 
plasma length scale L is approximately 25 cm, and the 
typical ion gyroradius is pi « 2 cm. The characteristic 
Alfven transit time is ta ~ 1-5 /is; plasma lifetime is 
about 15 /zs; resistive diffusion time is on the order of 
1 ms. The Lundquist number S = ^aLVx/ri ~ 10^-10"^, 
where Va is the Alfven speed and rj is the classical resis- 
tivity. Thus, the 15 /iS duration plasma dynamics lasts 
several Alfven times, and the magnetic flux is reasonably 
frozen into the plasma. 

III. EXPERIMENTAL RESULTS 
A. Plasma breakdown and initial evolution 

The Paschen condition for gas breakdown2SiS& shows 
that the voltage required for gas breakdown between two 
electrodes separated by distance d depends on the prod- 
uct pd, where p is the pressure in the region between 
electrodes. The dependence is such that the breakdown 
voltage becomes infinite when the pd product becomes 
less than about half the value at which the breakdown 
voltage is minimum. This property is exploited using a 
fast, localized gas injection, as shown in Fig. 3(bottom). 
The gas valve operation is adjusted so that pd satisfies 
the Paschen condition along the bias field lines. Because 
of the 6 mm size of the inter-electrode gap and the neg- 
ligible gas pressure in the gap, the pd value in the gap 
is so small as to be well to the left of the Paschen min- 
imum such that breakdown cannot occur there. This 
setup demonstrates the counter-intuitive fact that a very 
small distance between conductors in vacuum can provide 
perfect electrical insulation. This fact was also utilized 
in the Caltech solar prominence experimental setup. 

Breakdown along the bias field lines is shown in frame 1 
of Fig. 5, which is a CCD camera image of the plasma 
within approximately the first microsecond after break- 
down. The electrodes are located in the right hand side 
of each frame; the circular gap between electrodes is iden- 
tified by an arrow. The locations of the eight bright legs 
correspond to the locations of the eight pairs of gas in- 
jection holes, and the shape of the legs matches the cal- 
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culated bias field configuration. These observations indi- 
cate that breakdown is occurring along the bias field lines 
and that the gas pressure is non-uniform and concen- 
trated azimuthally at the eight positions corresponding 
to the gas injection holes. The breakdown occurs approx- 
imately 4-5 fis after high- voltage is applied (see Fig. 4). 
In order to minimize jitter and maximize reproducibility, 
the gas injection timing was adjusted to minimize the 
time delay before breakdown. 

The gun current ramps up steeply during the next 
10 fiS (see Fig. 4), and the eight bright legs expand and 
merge, as shown in the next two frames of Fig. 5. The 
rising gun current corresponds to an increase in toroidal 
magnetic flux linking the bias field lines, and therefore 
magnetic helicity is injected at the rate 2V^un?/'gun»2S. The 
applied gun voltage creates a radial electric field. From 
the radial component of the ideal Ohm's law, 

Ek + U^Bz - UzB^ = 0, (1) 

it is seen that a combination of toroidal rotation J/^ and 
axial flow Uz will arise to balance E'r. Note that an 
astrophysical accretion disk is similar to a coaxial gun 
but with the driving term coming from the Keplerian disk 
rotation C/^ which gives rise to i?R.^^ As described below, 
further evolution of the plasma after frame 3 of Fig. 5 is 
dependent on the peak value of Agun = Mo-^gun/V'gun- 



driven models of astrophysical jet coUimationjiSiiSiiSiii 

in which an ionized rotating accretion disk winds up a 
background magnetic fleld and injects magnetic helicity 
into the disk corona. 



2. Kink instability of central column at intermediate Agun 

For 40 m~^ < Agun ^ 60 m"-'^, the plasma column 
forms along the Z direction but then develops a helical 
instability with toroidal mode number n = 1, as shown 
in Fig. 8. The helical instability is shown to be consistent 
with an ideal MHD kink instability by two independent 
experimental measurements described below. 

The Kruskal-Shafranov condition for MHD kink insta- 
bility in cylindrical geometry is 

q{a) = 2TTaBz{a)/LB^{a) < 1, (2) 

where q is the safety factor and a and L are the column 
radius and length, respectively. By invoking the rela- 
tionships Tpgun ~ Tra^Bz{a) and B^{a) — /io-^gun/27ra = 
Agun^'gun/a, the Kruskal-Shafranov condition for kink in- 
stability can be expressed aai^/^ 



B. Plasma morphologies resulting from Agun 
parameter scan 

After the initial formation stage, the plasma can 
evolve into three distinct morphologies depending on 
peak Agun"^ Similar regimes were also observed in previ- 
ous Caltech experiments using a conventional cylindrical 
coaxial guniii^i The parameter Agun can be thought of as 
a boundary condition imposed on the plasma at the gun 
surface. Its functional form Agun = Mo^gun/V'gun can be 
easily derived by integrating V x B = AB over the inner 
gun electrode surface. The parameter Agun is varied ex- 
perimentally by adjusting the main capacitor bank volt- 
age (which controls /gun) and the bias coil bank voltage 
(which controls ipgun)- A plot of Agun parameter space 
is shown in Fig. 6, with the different observed plasma 
morphologies (I, II, III) indicated. 

1. Formation of stable plasma column at low Agun 

The formation of a central plasma column along the 
Z direction is observed for values of Agun ^ 40 m"^, 
as shown in Fig. 5. Magnetic probe measurements in 
this regime (Fig. 7) show that the column has magnetic 
field radial profiles resembling a screw pinch. The ax- 
ial expansion of the column is derived from Fig. 5 to 
be approximately 40 km/s which is of the order of the 
Alfven speed. The filamentary shape of this column pro- 
vides compelling experimental evidence for magnetically- 



Agun > 47r/L. (3) 

Thus, the experimentally observed instability onset can 
be compared to the Kruskal-Shafranov condition by 
knowing only Agun and L. Figure 9 shows a scatter plot 
of Agun versus L for an ensemble of shots in which a 
plasma column forms. The data points represent both 
diflerent shots and different times within the same shot 
{i.e., before and after the instability occurs). It is seen 
that the threshold for the appearance of the helical insta- 
bility is in good agreement with the Kruskal-Shafranov 
limit (dashed line in Fig. 9). 

Consistency with the Kruskal-Shafranov limit was also 
checked with direct magnetic probe measurements. Ra- 
dial proflles of Bz and B^ yielded q profiles as a function 
of time, as shown in Fig. 10. It is shown that q drops to 
unity around 10.5 /xs. The corresponding image sequence 
for this shot, shown in Fig. 11, indicates that the helical 
instability develops at this same time. Thus the observed 
helical instability is consistent with an ideal kink, with 
onset occurring when q drops to unity. The q is initially 
greater than unity because the bias poloidal fiux domi- 
nates and L is small. Then L increases quickly and B^j, 
also increases to a lesser degree, and therefore q drops 
according to Eq. 0. 

The growth of the kink mode was captured by short- 
ening the inter-frame time of the CCD camera to 0.25 /xs 
(see Fig. 12). The mode amplitude as a function of time 
is shown in Fig. 13 and is seen to exhibit a fairly linear 
growth rate. 
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3. Quick plasma detachment at high Agun 

At values of Agun 60 ni^^, the plasma appears to 
"pinch-ofT" without the formation of a central column, 
as shown in Fig. 14. From the image sequence, the prop- 
agation velocity of the detached plasma along the Z-axis 
is calculated to be around 50 km/s (of the order of the 
Alfven speed). Magnetic probe measurements indicate 
that the detached plasma has much stronger toroidal 
than poloidal field (Fig. 15). It is possible that this 
toroidal-flux-rich detached plasma will relax to a sphero- 
mak configuration with equivalent toroidal and poloidal 
flux. 

Poloidal flux amplification, as measured by the mag- 
netic probe array, is large in this regime. Flux amplifica- 
tion for all three regimes is shown in Fig. 16 as a function 
of peak Agun- In the kinked regime, the flux amplification 
is typically around a factor of two, and the flux amplifica- 
tion mechanism in this regime (discussed below) has been 
identified. The flux amplification is neligible in the low 
Agun regime, but it is very large in the high Agun regime. 
The mechanism for the high Agun regime is not known, 
but it is probably related to the mechanism operating 
in the kink regime. An interesting question requiring 
further study is what determintes the maximum amount 
of flux amplification. The next section shows that the 
kink is a mechanism for poloidal flux amplification and 
spheromak formation. 



tion, i.e., k • x = constant, such that 

A:^^ + fc^i?!/) = constant. (5) 

Thus, substituting Eq. Q into Eq. ISJ gives the coordi- 
nates of the helix, (j> = {B^/ RBz)Z + constant, from 
which it can be seen that the kinked current channel 
is a right-handed helix if JzBz ^ B^Bz > and a 
left-handed helix if JzBz < 0. Thus, the helix will 
always have the form of a solenoid with the appropri- 
ate handedness to amplify the original Bz. The addi- 
tional ip introduced by the kink can be estimated by con- 
sidering the helically deformed current channel to be a 
solenoid with current /, turns per length 1/L, and radius 
a. The solenoid formula for the field inside the solenoid 
is Bz — fioI/L; the produced by the solenoid is ira^Bz 
and thus depends non-linearly on the kink amplitude a. 
Using the measured values a « 5 cm, L w 20 cm, and 
/ « 60 kA at 13.5 /is, the -0 generated by the kink is 
predicted to be approximately 1 mWb, which is within 
a factor of 2 of the observed amplification of V'max over 
V'gun- The discrepancy is within the accuracy of a and 
L measurements and of the ip calculation assuming ax- 
isymmetry in the presence of the rotating kink. Because 
the coaxial gun injects only toroidal flux, 3D plasma dy- 
namics must be responsible for ijj exceeding '0gun- The 
dynamics are provided by the kink, which is the mecha- 
nism by which toroidal flux is converted to poloidal flux. 



C. Kink as mechanism for poloidal flux 
amplification and spheromak formation 

The kink instability is identified experimentally as 
a mechanism for (1) poloidal flux amplification and 
(2) spheromak formation. Both are experimentally ob- 
served immediately following kinking of the central col- 
umn. This can be seen by examining the evolution of 
poloidal flux contours and maximum ■0 (Fig. 17) calcu- 
lated from time-resolved magnetic probe data and the 
sequence of CCD images showing the kink (Fig. 11). 

1. Paramagnetism of the kink 

The kink modifies the direction of current flow from 
purely Z (poloidal) to a partially cf) (toroidal) direction. 
Consequently, it converts toroidal to poloidal flux. This 
process is paramagnetic, amplifying ip over the initial 
applied i/'gun = 1-7 mWb. The paramagnetism is un- 
derstood by realizing that kinks involve perturbations 
with dependence exp(zk ■ x) where k • B = 0. The latter 
means 

fcz = -k^B^/Bz, (4) 

where A:r = 0. The trajectory of the kink is determined 
by considering a locus of constant phase of the perturba- 



2. Evidence for spheromak formation 

The appearance of the kink is followed immediately by 
three signatures of spheromak formation: (1) appearance 
of closed -f/j-contours [calculated assuming axisymmetry, 
^{R,t) = J^2nR'Bz{R',t)dR'], (2) V-amplification, 
and (3) magnetic field radial profiles consistent with 
spheromak formation. It is important to note that the 
relationship between closed ^/'-contours and closed flux 
surfaces becomes ambiguous when axisymmetry is bro- 
ken. Thus, in the presence of a non-axisymmetric kink, 
closed T/i-contours indicate closed flux surfaces only in 
a time-averaged manner. Shortly after 13 /xs, the kink 
breaks apart (Fig. 11). This coincides with signatures of 
spheromak formation as observed in the magnetic probe 
measurements (Fig. 17). At approximately 13 /iS, closed 
i/i-contours appear and V'max is amplified to larger than 
"Vgun ("V-a-mplification is due mainly to broadening of the 
Bz profile after 12 /is). At 15 /iS, magnetic field pro- 
files consistent with spheromak formation are observed 
(Fig. 18). The measured radial profiles of Bz and B^ at 
15 /iS are compared with Taylor state solutions^ in cylin- 
drical geometry, i.e. uniform A solutions of V x B = AB. 
The solutions are Bz ~ Jo{XR) and B^ ^ Ji{XR), where 
Jo and Ji are Bessel functions of order zero and one, 
respectively, and the best fit is found for A « 15 m~^ 
with a radial offset of 4 cm (displacement of spheromak 
off the geometric axis). The slight disagreement between 
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measured profiles and the Taylor solution is not surpris- 
ing since the spheromak is expected to be either (1) still 
undergoing relaxation toward the Taylor state or (2) in 
a modified relaxed state since it is still being driven by 
the gun, where peak Agun ~ 50 m~^. 



IV. DISCUSSION 
A. Spheromak formation and topology 

Spheromak formation involves a sequence of steps that 
in principle should be possible to visualize mentally. 
However, the inherent geometric and topological com- 
plexity of these steps has made such a visualization so 
challenging that no consensus exists on the precise form 
of this sequence. This shortcoming is important be- 
cause without being able to visualize how the 3D topol- 
ogy evolves during the spheromak formation sequence, 
it is difficult to optimize this process and the closely re- 
lated sustainment process. "^^ Because visualization of the 
process has proved difficult, the main appoach used to 
date for experimentally characterizing spheromak forma- 
tion and sustainment is modal analysis, a method where 
Fourier modes are first identified and then the time de- 
pendence and spatial profile of these modes are measured 
and compared with theoretical models. While useful in 
many respects, modal analysis gives little insight into 3D 
topological issues such as linkage, connectivity, knotted- 
ness, and conversion between toroidal flux and poloidal 
flux. The photographs of evolving magnetized plasmas in 
the present experiments have provided an actual visual- 
ization which helps provide new insights into these issues. 
New interpretations are given here and placed into con- 
text with respect to past experimental workiSiSSi^iiSSiSi^i 
It should also be noted that many aspects of the present 
work are qualitatively similar to recent time-dependent 
MHD simulations of spheromak formation via electro- 
static helicity injection. 

Bounded, axisymmetric, force- free plasma equilibria 
were first discussed in the 1950's.^i^S*^£i2& These theo- 
retical models invoked rather abstract astrophysical con- 
texts and did not contain descriptions of the detailed 
processes by which such equilibria could be created. A 
hypothetical analogy to this situation would be to be in 
possesion of an abstract theoretical model demonstrating 
that soap bubbles should in principle exist, while never 
having observed a soap bubble nor being aware of how 
to blow soap bubbles. The soap bubble analogy is apt 
because a bounded, axisymmetric, force-free equilibrium 
can be thought of as being a bubble (local maximum) of 
poloidal flux iIj{R,Z). First derivatives (gradients) of 
show that a local ip{R, Z) maximum constitutes a mag- 
netic axis, and second derivatives (Laplacian-like opera- 
tors) of ip show that a toroidal current necessarily flows 
along this magnetic axis. 

Lacking detailed knowledge for how to make these 
"magnetic bubbles" is a non-trivial shortcoming because 



Cowling's anti-dynamo theorem^ shows that there is no 
simple way to make such bubbles. Specifically, Cowl- 
ing showed that creation/sustainment of an axisymmet- 
ric, steady-state toroidal current cannot be achieved via 
axisymmetric {i.e., simple) means. Since creation of an 
axisymmetric toroidal current is tantamount to creating 
a local ip{R, Z) maximum, the Cowling anti-dynamo the- 
orem effectively requires that any poloidal flux amplifl- 
cation process must be non-axisymmetric. Thus, sphero- 
mak formation and sustainment necessarily involve non- 
axisymmetric processes. 

Helical {i.e., non-axisymmetric) distortions of magne- 
tized coaxial gun central columns and associated poloidal 
flux amplification were first observed in the early 1960's 
in experiments by Alfven et ali^ and Lindberg et ali^ 
It is interesting to note that these early magnetized 
coaxial gun experiments revealed phenomena relating to 
spheromak formation nearly two decades before the word 
"spheromak" was coined by Rosenbluth and Bussac4£ In 
particular, Lindberg et al»^ imaged the helical instability 
of an initially azimuthally symmetric axial current chan- 
nel, characterized this as an "m = 1 instability," noted 
that "such instabilities are well known in the theories of 
pinches," and measured an amplification of the poloidal 
flux. Lindberg et al. interpreted this central column he- 
lix as a kink instability of the geometric axis current and 
postulated that the helical deformation of the current 
path was responsible for observed poloidal flux amplifica- 
tion. However, Lindberg et al. did not provide quantita- 
tive measurements establishing that the observed helical 
instability was indeed a kink {i.e., that the configuration 
destabilized when the q = I kink stability threshold was 
crossed) nor that the instability produced a magnitude 
and polarity of poloidal flux consistent with observation. 

Similar modes (now called n ~ 1 modes) were also 
observed in all later coaxial gun experimentsri*224Siii 
but there have been two distinct, differing interpretations 
of these modes. In the first interpretation, referred to 
here as "rotating static mode," the modes are assumed 
to have a steady-state geometry in some rest frame so 
that their apparent laboratory frame time-dependence is 
solely due to the rotation of the mode rest frame with re- 
spect to the laboratory frame. In this interpretation the 
modes are thus assumed to have a laboratory frame time- 
dependence ^ cos(n6' — cot) so that an observer located 
in a frame rotating with the mode {i.e., a frame with 
coordinate 9' — 9 — wt/n) would see a time- independent, 
or static, mode. In the second interpretation, referred to 
here as "relaxation oscillation," the modes are not static 
in any frame but, rather, are growing modes which start 
with zero amplitude, build up to finite amplitude, un- 
dergo a nonlinear crash, and then the process repeats. 
There might be rotation as well, but even if one moved 
to an optimally chosen rotating frame, there still would 
be relaxation oscillations. Thus, in the relaxation oscil- 
lation case, the mode time-dependence during its growth 
stage would be ~ e"*"* cos n9 or possibly ~ e"*"* cos{n9 — u!t) 
so that no change of observer rest frame would cause the 
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mode to appear static. 

The n = 1 modes were typically observed as oscilla- 
tions in magnetic field during the sustainment (as op- 
posed to formation) phase. However, it seems likely that 
sustainment ought to involve much the same processes 
as formation since formation requires creating a localized 
poloidal flux maximum, and sustainment requires main- 
taining this maximum against resistive losses. Knox et 
al^ observed n — 1 and higher order modes using wall 
magnetic probe arrays and proposed that sustainment 
was due to a rotating static mode which they called a 
"rotating internal kink distortion." On the other hand 
Nagata et al^ made similar magnetic probe measure- 
ments but interpreted sustainment as a sequence of relax- 
ation oscillations ( "collapse and recovery" ) . Thus Knox 
et al. considered the mode as a rotating perturbation of 
an otherwise axisymmetric flux surface, whereas Nagata 
et al. interpreted the mode as being the creation and de- 
struction of flux surfaces. 

The SPHEX groupS2i22i^ used internal probes to iden- 
tify two distinct structures in their experiment, namely 
(i) a central column of open magnetic field lines with a 
force-free current flowing from the gun electrode along 
the geometric axis and (ii) an annulus (torus) of closed 
field lines and toroidal current linking the central column. 
The SPHEX group found that the n — 1 oscillations 
were primarily in the central column and, using probes, 
showed that the central column had a helical deforma- 
tion. The helical deformation was inferred on the as- 
sumption that the n = 1 oscillation was a rotating static 
mode, and this mode was interpreted to be a nonlinearly 
saturated kink instability of the geometric axis current 
rotating toroidally at the E x B frequency.''^ It is impor- 
tant to emphasize that the rotation was assumed^^ but 
not directly measured. Because the signal was observed 
to be periodic, it was assumed that the fluctuations were 
of the form cos{d—ujt) so that observations made at differ- 
ent times could be used to map out the spatial proflle of 
the assumed rotating structure. However, this assump- 
tion of rotation misses the interpretation that there is 
a relaxation oscillation with repeated births of new he- 
lices, i.e., the rotation assumption did not consider the 
possibility of a repetitive sequence of modes with time 
behavior e''* cos 6. 

The present results suggest that the n = 1 fluctu- 
ations result from repeated births of new kinked he- 
lices with time dependence ^ e^'cos^ and that these 
helices detach and merge into the time-averaged closed 
flux region of a spheromak. According to this interpre- 
tation of spheromak formation/sustainment, the current 
starts out aligned along the geometric axis and so is ini- 
tially purely poloidal. Then, upon kinking (i.e., the cur- 
rent channel becoming helical so that it is effectively a 
solenoid), the current channel axis has a toroidal com- 
ponent which creates new poloidal flux. Kink growth 
means that the radius of the helix increases. Growth 
of the kink, in turn, means that the radius of the he- 
lix increases until this radius becomes so large that the 



helical current channel merges with any existing sphero- 
mak toroidal current, thereby amplifying/sustaining the 
poloidal flux. New current will now flow from the gun 
along a new straight central column aligned along the 
geometric axis, and initially, this new central current col- 
umn will be stable against kinking because of the strong 
axial magnetic field resulting from the strong toroidal 
current flowing in the spheromak. However, because no 
electromotive force sustains the spheromak toroidal cur- 
rent against Ohmic decay, the spheromak toroidal cur- 
rent decays with time so that the poloidal field on the 
geometric axis becomes weaker. The open poloidal field 
on the geometric axis weakens to the point that a new 
kink develops on the geometric axis. The new kink de- 
forms into a new solenoid which merges with the sphero- 
mak, replenishing the toroidal curent in the spheromak, 
and increasing the geometric axis open poloidal field so 
that the process starts over again. Thus, power flow from 
the geometric axis region to the spheromak is not the re- 
sult of rotation of a saturated kink but is instead due to 
the repetitive creation of successive helical (kinked) ge- 
ometric axis current channels which expand nonlinearly 
in radius and merge with the spheromak. 

The interpretation in the above paragraph is but- 
tressed by considering what happens when a simple 
solenoid, consisting of many turns of wire, is rotated 
about its own axis. Such a rotation does not result in 
flow of magnetic energy or in transformation of toroidal 
field energy into poloidal field energy. In fact rotation of 
a solenoid about its axis does not change the magnetic 
field produced by the solenoid at all because rotation of 
a solenoid does not change the relative toroidal drift ve- 
locity between electrons and ions in the solenoid wire. In 
contrast, increasing the radius of the solenoid, as would 
be the consequence of a growing kink instability, increases 
the poloidal fiux through the bore of the solenoid. If the 
conductor forming the solenoid {i.e., the initial geomet- 
ric axis current channel) is a flux conserver, then this 
increase of poloidal flux linking the conductor must be 
at the expense of a decrease in the toroidal flux linking 
the conductor, resulting in a conversion from toroidal to 
poloidal flux due to the kink. 

The primary result in this work regarding spheromak 
formation is the experimental demonstration that the 
kink instability of the geometric axis current and its 
growth result in the poloidal flux ampliflcation required 
to produce the local maximum of poloidal flux necessary 
for spheromak formation. This result is qualitatively sim- 
ilar to the point of view expressed by Lindberg et al. and 
Nagata et al. but differs signiflcantly from the point of 
view expressed by Knox et al. and the SPHEX group. 
The present work has provided a quantitative demonstra- 
tion that the onset of the helical instability is consistent 
with the Kruskal-Shafranov condition for kink instability 
and a theoretical demonstration that the kink is necessar- 
ily paramagnetic {i.e., amplifies poloidal fiux). Further- 
more, it was shown that the poloidal flux ampliflcation 
calculated using the visually observed kink dimensions 
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are consistent with the magnetic probe measurements of 
poloidal flux. These measurements support the inter- 
pretation that poloidal flux amplification is due to kink 
growth and not kink rotation. Repetition of this sequence 
with successive, new kinked current channels during sus- 
tainment could explain the previously observed n = 1 
fluctuations. The flux amplification mechanism would be 
insensitive to the magnitude or direction of the rotation 
frequency. 



B. Astrophysical jets 

Magnetically-driven astrophysical jets are intriguing 
astronomical phenomena which have been observed for 
nearly a century. These jets are highly coUimated and 
have been observed to emanate from stars, black holes, 
and active galactic nuclei; the sources are generically 
called central objects. The characteristic size of jets 
ranges over many orders of magnitude, but in all cases 
the axial extent of the jet is much larger than the charac- 
teristic dimension of the source. Theoretical models of jet 
formation and collimation have been discussed at length 
in the literatureiiS'ii'i^ Observations of helical magnetic 
fields in jeta^ provide especially tantalizing evidence that 
plasma physics plays an important role governing jet col- 
limation, structure, and dynamics. Theoretical models 
based on MHD have become the leading candidates to 
explain jet formation and collimation. Many of these 
models show how axisymmetric force- free or nearly force- 
free plasma equilibria could be made to resemble the 
long thin morphology of a iet^i^^^ More recently, time- 
dependent 3D MHD simulationsiSii^ have illustrated dy- 
namical collimation starting with an initial "seed" mag- 
netic field configuration. Despite impressive similarities 
between the theoretical models and astronomical obser- 
vations of jets, experiments to test the underpinnings of 
these theories have been lacking. 

The planar spheromak gun formation experiment dis- 
cussed here has a geometry very similar to the geome- 
try associated with the formation region of magnetically- 
driven astrophysical jets (Fig. 2), and so this work can be 
considered as a laboratory investigation into the underly- 
ing processes governing astrophysical jet formation, colli- 
mation, and accelerationii4 The MHD theories are based 
on the generation of an axisymmetric poloidal current 
with geometry essentially the same as the poloidal cur- 
rent in the central column of the present experiments. In 
the astrophysical situation, this poloidal current is driven 
by a radial electric field resulting from the rotation of 
an accretion disk located in the equatorial plane of the 
central object. The rotation of the accretion disk cuts 
across a pre-existing poloidal magnetic field, resulting in 
an electric field E = —U^cj) ^ BzZ = —UtpBzR, which is 
equivalent to the radial electric field applied across the 
gap between inner and outer electrodes in the experi- 
ment. In both cases, magnetic helicity is injected at the 
rate dK/dt = 2Vip, where V is the potential drop be- 



tween the inner and outer electrodes (or between star 
and disk) and is the poloidal magnetic flux linking the 
two. 

The present work demonstrates experimentally for the 
first time that magnetic helicity injection in relevant ge- 
ometry leads to the formation and collimation of a jet 
(Fig. 5). In particular, the data indicate both a jet-like 
fiow and a collimation of this flow, i.e., the lengthening 
central column is the jet, and the narrowing of the cen- 
tral column is the collimation. The jet flow is consistent 
with a recent model^^ proposing that axial gradients in 
drive axial flows from where the current channel ra- 
dius is small to where the channel radius is large and that 
collimation occurs where V • U is negative so that there 
is a compression of the toroidal flux frozen into the flow. 
The present work also shows that collimation can occur 
within just one twist of the poloidal magnetic field, as 
demonstrated by the fact that the kink instability occurs 
when q drops to unity, after collimation has already been 
observed. It is interesting to note that the kink insta- 
bility has been proposed to explain the observed knotted 
structures in astrophysical jets4& 



V. SUMMARY 



Magnetic helicity injection experiments using a simple 
magnetized coaxial gun setup have revealed insights into 
the dynamics of spheromak formation and collimation of 
magnetically-driven astrophysical jets. Depending on the 
peak value of Agun, three different plasma morphologies 
were observed: (1) a collimated jet at Agun ^ 40 m~^, 
(2) kink instability of the jet at 40 < Agun S 60 m~^, 
and (3) a quickly detached plasma at Agun ^ 60 m~^. 
The kink instability was shown to provide the mechanism 
for poloidal flux amplification in coaxial gun spheromak 
formation. The latter result supports a model in which 
repeated births of kinked helices can form and sustain 
spheromaks. It is proposed that this mechanism, rather 
than a rotating static saturated kink which does not lead 
to flux ampliflcation, could explain past observations of 
n = 1 fluctuations in coaxial gun spheromak experiments. 
A collimated jet is formed by driving electric current 
along a background poloidal held in a geometry relevant 
for magnetically-driven astrophysical jets. Collimation is 
observed well before kink onset. Since kink onset takes 
place when the magnetic field has undergone one com- 
plete twist, consistent with the Kruskal-Shafranov the- 
ory, these observations demonstrate that collimation oc- 
curs before the initial poloidal field has undergone one 
complete twist. Kink dynamics could be related to wig- 
gled structures in astronomical observations of galactic 
jets. 
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Figure captions 

Figure 1. End-view photograph of coaxial gun and 
rotatable magnetic probe array (end-view). 

Figure 2. Schematic of essential ingredients of both a 
coaxial gun and astrophysical star-disk system. 

Figure 3. Top: large vacuum chamber with multiple 
diagnostic ports and coaxial gun mounted on right end- 
dome. Bottom: Side-view schematic of planar coaxial 
gun, showing inner electrode (blue), outer electrode 
(green), gas feed lines, contours of constant bias poloidal 
flux, and cylindrical coordinate system. 

Figure 4. Typical /gun (solid) V^un (dashed) traces. 
Breakdown occurs a.t t = 2-3 fxs when V^un drops and 

-/gun riSGS • 

Figure 5. CCD images of plasma evolution (shot 1210) 

for low Xgun regime, in which a central plasma colmim 
forms. The circular gap between electrodes is indicated 
by an arrow. 



Figure 16. Flux amplification versus peak Agun- 

Figure 17. Top: poloidal magnetic field vectors and ip 
contours as a function of R and t. Bottom: Time evolu- 
tion of at i? = 19 cm (shot 2472). Flux amplification 
is observed during time of the kink and appearance of 
closed tp contours. 

Figure 18. Comparison of Bz (diamonds) and 
(squares) profiles at i = 15 /zs with Taylor state (solid 
line). Profile location is indicated by vertical dashed line 
in top panel of Fig. 17. 



Figure 6. Plot of Agun = Mo-^gun/V'gun parameter space 
showing the three plasma regimes. 



Figure 7. Magnetic field and total current as a function 
of R for low Agun regime central column (shot 2465). 

Figure 8. CCD images of plasma evolution (shot 1233) 
for intermediate Agun kinked plasma regime. 

Figure 9. Plot of Agun versus column length L, 
showing consistency with Kruskal-Shafranov condition 
Agun > 47r/L for MHD kink instability. 

Figure 10. Radial profiles of g = 2'kRBz/LB^ (shot 
2472). The column develops a helical instability at the 
same time as q near the axis drops to unity, also showing 
consistency with Kruskal-Shafranov condition for MHD 
kink instability. 

Figure 11. Images of kink development (shot 2472). 

Figure 12. High time-resolution images of kink growth 
(shot 1247). Inter-frame time is 250 ns. 



Figure 13. Mode amplitude of kink perturbation versus 
time. 



Figure 14. Images of plasma evolution at high A, 
(shot 1181). 

Figure 15. Radial profiles of Stor and Bp^i for high A, 
plasma at t = 12 fis (shot 2457). 
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